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Abstract 

One  way  to  perform  CAA  simulations  is  to  split  the  flow  field  in  a steady  mean  flow  and  turbulent 
and  acoustic  fluctuations.  The  small  acoustic  fluctuations  can  then  be  calculated  from  linearized 
Euler  equations  (LEE).  The  mean  flow  can  be  supplied  by  the  solution  of  the  Reynolds  averaged 
Navier-Stokes  (RANS)  equations.  Since  different  numerical  requirements  exist  for  the  numerical 
grids  of  RANS  and  LEE  simulations,  an  interpolation  procedure  between  RANS  and  LEE  grids  is 
necessary.  The  following  paper  presents  an  interpolation  procedure  between  different  structured 
multi-block  grids  and  its  application  to  RANS/CAA  simulations  of  a generic  wing  section. 


1.  Introduction 

The  most  accurate  way  to  calculate  the  generation  and  propagation  of  sound  in  flows  is  the  solution  of  the  unsteady 
conservation  laws  of  mass,  momentum,  and  energy.  Since  this  direct  numerical  simulation  is  even  for  simple  flows 
an  extremely  expensive  task,  a lot  of  approximative  methods  have  been  proposed  in  order  to  attack  real  world 
problems.  One  of  them  is  the  splitting  of  the  flow  field  in  mean  flow  and  turbulent  and  acoustic  fluctuations.  The 
mean  flow  is  approximated  by  a solution  of  the  Reynolds  averaged  equations  (RANS),  the  turbulent  fluctuations 
e.g.  by  an  SNGR  model  (Stochastic  Noise  Generation  and  Radiation)  [2,3,4,7],  and  the  acoustic  fluctuations  by  the 
solution  of  the  linearized  Euler  equations  (LEE).  This  approach  is  utilized  in  the  german  SWING  project 
(Simulation  of  Wing-flow  Noise  Generation),  which  shall  provide  a general  computer  code  for  the  calculation  of 
airframe  noise,  especially  that  generated  by  high  lift  devices. 

Using  this  RANS-LEE  approach,  it  is  generally  necessary  to  use  different  grids  for  both  calculations.  The  reason 
for  this  are  the  different  numerical  requirements  in  both  cases.  E.g.  the  RANS  calculations  need  very  fine  grids  in 
the  boundary  layers  of  the  body,  whereas  the  LEE  require  a sufficient  fine  resolution  in  the  farfield  of  the  body. 
Therefore  the  RANS  solution  has  to  be  provided  on  the  CAA  grid  by  means  of  interpolation.  The  present  paper 
discusses  this  transfer  of  RANS  data  into  the  CAA  perturbation  code  in  view  of  robustness,  accuracy  as  well  as  the 
numerical  implications  on  CAA  simulations  of  airframe  noise. 

This  article  is  structured  as  follows.  First  a brief  overview  of  RANS  calculations  of  a generic  wing  section  will  be 
given.  Then  the  CAA  grid  and  the  interpolation  procedure  between  the  RANS  grid  and  the  CAA  grid  will  be 
considered.  Finally  some  results  of  CAA  calculations  with  a Navier-Stokes  and  Euler  mean  flow  will  be  presented. 

2.  The  RANS  Calculations 

As  generic  example,  the  viscous  flow  around  a wing  section  of  an  airfoil  representative  for  those  at  modem  civil 
aircraft  was  calculated.  The  flow  solver  was  the  DLR  FLOWer  code  [1,8,9],  which  solves  the  compressible 
Navier-Stokes  equations  oh  a grid  consisting  of  topologically  rectangular  blocks  (structured  multi-block  grid).  The 
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free  stream  Mach  number  was  Ma  = 0.2  , the  Reynolds  number  Re  = 1 .6  x 1 06  and  angle  of  attack  of  (X  = 7° . As 
turbulence  model,  the  k- to  model  of  Wilcox  was  used. 


2.1  The  RANS  Grid 

The  RANS  grid  consists  of  12  blocks  with  a total  of  about  86000  grid  points,  cf.  Fig.l.  The  topologically 
complicated  structure  of  the  grid  stems  from  the  fact  that  with  the  same  topology  also  a wing  section  with  deployed 
Fowler  flap  has  to  be  calculated  (as  part  of  a future  3d  simulation).  It  can  be  seen  that  a typical  RANS  grid 
becomes  very  coarse  in  the  far  field  region.  Any  calculation  of  wave  propagation  would  break  down  there. 
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Fig.  1 The  RANS  grid(detail).  (Thick  lines  are  block  boundaries) 

2.2  The  Mean  Flow  Field 

The  mean  flow  field  around  a wing  section  at  high  angle  of  attack  has  a strong  pressure  minimum  at  the  nose  of 
the  profile,  cf.  Fig.3,  which  corresponds  to  high  local  velocities  and  gradients  of  the  flow  variables.  The  maximum 

Mach  number  at  the  nose  is  about  Ma  = 0.45  in  contrast  to  a free  stream  mach  number  of  0.2 . The  thin 
boundary  layer  requests  a small  grid  spacing  there.  In  the  present  case  the  distance  of  the  first  grid  line  from  the 

surface  was  1 x 1 0 5 of  the  chord  length.  The  examined  wing  section  has  a blunt  trailing  edge  with  a thickness  of 
about  0.33%  of  die  chord  length.  This  finite  trailing  edge  requests  a considerable  amount  of  fine  tuning  of  the  grid 
in  order  to  achieve  sufficient  convergence  of  the  RANS  solution.  Besides  problems  with  the  turbulence  models 
there,  the  main  reason  for  this  difficulty  is  the  complex  grid  topology  in  the  wake  of  the  profile.  For  comparison 
purposes  also  an  Euler  calculation  was  performed  which  needed  a different  grid  not  shown  here. 
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3.  Interpolation  between  RANS-  and  CAA- 
Grid 


3.1  The  CAA  Grid 

The  CAA  grid  in  the  present  example  consists  of  2 
blocks  with  a total  of  about  1 8000  grid  points,  cf.  Fig.3. 
The  boundary  layers  of  the  profile  were  resolved  by 
approximately  10  grid  points  with  a minimum  grid 

spacing  of  5 xl0~6  at  the  nose  of  the  profile.  This 
results  in  an  allowable  time  step  of  the  CAA  calculation 
of  the  same  magnitude.  The  minimum  grid  spacing  at 
the  nose  is  about  the  half  of  the  spacing  of  the  RANS 
grid  there.  This  small  spacing  was  chosen  mainly  in 
order  to  resolve  possible  small  unsteady  flow  structures. 
It  should  be  remarked  however  that  the  progression  of 

Fig.  2 The  pressure  distribution  of  the  mean  flow  **  31  the  nose  is  larSer  ^ the  progression 

choosen  in  the  RANS  grid.  The  gridspacing  in  the  far 
field  was  chosen  such  that  one  retains  enough  resolution  of  the  evanescent  sound  waves  and  is  therefore 
significantly  smaller  than  in  the  RANS  grid. 


Fig.  3 The  CAA  grid  (2  blocks.  Thick  lines  are  block  boundaries.) 

It  should  be  noted  as  well,  that  a considerable  coarsening  (factor  20)  of  the  CAA-grid  normal  to  the  boundary  layer 
close  to  the  surface  will  still  allow  for  an  adequate  representation  of  essential  features  of  the  physics.  As  will  be 
seen  later,  this  is  a consequence  of  the  fact  that  the  presence  of  boundary  layers  has  only  a small  influence  on  the 
scattering  problem  considered. 
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3.2  The  Interpolation  Procedure 

The  interpolation  procedure  between  the  RANS  grid  and  the  CAA  grid  is  complicated  by  two  circumstances.  First, 
the  two  grids  can  be  located  almost  arbitrary  in  3d  space  (The  only  restriction  is,  that  the  CAA  grid  has  to  be 
enclosed  by  the  RANS  grid).  This  results  in  an  expensive  search  for  appropriate  initial  conditions  for  the  iterative 
determination  of  the  parameters  of  a point  of  the  CAA  grid.  The  second  issue  is  the  very  fact  that  in  complicated  3d 
calculations,  one  has  to  search  some  million  CAA  points  in  some  million  RANS  cells.  This  poses  great  demands 
upon  the  speed  of  the  algorithm.  Therefore  an  easy  to  evaluate,  local  polynomial  approximation  was  chosen  for  the 
representation  of  the  grid  functions. 

3.2.1  The  Interpolating  Polynomial 

The  RANS  solver  FLOWer  developed  at  DLR  utilizes  a general  curvilinear  multi-block  structure.  Every  block  of 
the  grid  consists  of  the  n,x«2x  «3  points 


, 0 </<«,- 1 , 0 < j < n2  - 1 , 0 < k < «3  - 1 


Every  flow  variable  in  the  FLOWer  code  is  defined  at  these  grid  nodes.  It  is  convenient  to  consider  the  grid 
coordinates  and  every  other  function  defined  on  the  block  as  function  of  the  grid  indices.  I.e.,  we  introduce 
parameters  (£,T|,£),  in  such  a way  that  e.g. 


x(i,j,k)  = 


x>jk 

y.jl: 

l Z„L 


Mi,  j,  k)  = x{^  = i,J]  = j,  £ = k) 

In  the  following,  the  symbol  / stands  for  any  variable  that  is  defined  at  the  grid  nodes.The  interpolation  problem 
can  now  be  formulated  as  follows 

1 . Choose  an  appropriate  representation  of  a function  f (£,  7],  £)  on  the  block. 

2.  Determine  for  a given  point  y inside  the  block  the  values  (£,  7 J,  from  the  vector  equation  y = x(£,  TJ,  £)  . 

Since  one  has  to  interpolate  a vast  amount  of  points  (up  to  some  millions  in  3d  calculations),  one  needs  a fast 
interpolation  method  with  acceptable  accuracy.  This  leads  to  a polynomial  interpolation  of  the  function  f on  every 
cell  of  the  block.  Therefore,  we  consider  now  a cell  i<£  <i  + 1 , j<TJ<  j + 1 , k < £ < k + 1 of  the  block  and 
introduce  local  coordinates  u=£  — i,v  = T]  — j,w  = £ — k in  this  cell.  The  interpolation  problem  now  boils 
down  to  a polynomial  interpolation  of  a function  on  a unit  cell  0<w<l,0<v<l,0<w'<l,  where  the  function 
values  (and  perhaps  appropriate  derivates  of  the  function)  are  given  at  the  comers  of  this  unit  cell.  The  most  simple 
interpolation  is  the  linear  one 

f(u,v,w)  = am+al00u  + a0l0v  + amw+au0uv  + at0luw+amvw+Gluuvw 

The  eight  coefficients  am, ,am  can  be  determined  by  the  eight  function  values 

f 000  = /(0,0,0), ,/ni  =/(l,l,l)  at  the  comers  of  the  cell.  Unfortunately  the  accuracy  of  this  interpolation  is 

quite  poor.  An  improvement  however  can  be  achieved  using  a higher  order  polynomial.  In  this  case  one  needs 
additional  information  from  the  neighboring  cells  or  from  derivatives  at  the  comers  of  the  cell  in  order  to  determine 
the  higher  number  of  coefficients.  In  order  to  construct  a local  approximation  the  second  approach  was  choosen. 


USing  aPProPnate  values  for  the  2 nd  ■ 

to  third  order  in  u , v and  *,  as  merpolatgZL(n  ^ ^ C°nVenient  consider  the  Mow 
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lng  polynomial  up 


f(»,v,w)~a 


000  +ai,,uvw+  agii 

U 


Now  one  has  to  detennine  the  32  coefficients  V—  y 

2a2cx>+6a3mU  + 2a  . o 

*“■  *■ -e  coefficie™,  + &™“*+K,,v„+6<Im„w 

(Zero  entries  are  mated  win,  ( )>  ^vely  fe,m  fc  f„,lowin 

(2  ■ . Va  s S recursively  solvable  linear  systems 
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3.22  Calculation  of  the  Parameters 

Once  an  interpolating  polynomial  can  be  constructed  for  any  function  f(£,T],£) , the  main  problem  of  the  grid  to 
grid  interpolation  can  be  tackled,  i.e.  the  determination  of  the  parameter  values  (£,  t],  for  an  arbitrary  point  y 
inside  the  RANS  grid.  This  can  be  split  into  two  parts 

1.  Find  the  appropriate  block  of  the  RANS  grid  and  some  initial  conditions  for  the  parameter  values. 

2.  Solve  the  (nonlinear)  vector  equation  y = x(£,  7],  £)  for  (£,  7],£)  using  the  polynomial  approximation  above 
for  every  component  function  of  x . 

The  crucial  point  in  this  algorithm  is  to  provide  appropriate  initial  conditions  for  the  parameters.  First,  y is  checked 
against  the  bounding  box  of  a block.  Ify  is  inside  the  bounding  box,  initial  conditions  are  searched  in  every  cell  of 
the  block.  In  order  to  do  this,x(^,r/,^)  is  approximated  on  every  cell  by  a linear  lunction.  Initial  values  for 
can  than  be  found  by  solution  of  a simple  system  of  linear  equations.  This  procedure  works  for  most  of 
the  cells  of  the  RANS  grid.  If  however  the  faces  of  a cell  have  strong  curvature  and/or  if  the  the  aspect  ratio  of  the 
cells  takes  very  large  values  this  procedure  can  fail.  This  can  happen  especially  with  cells  in  the  boundary  layer  at 
the  nose  of  the  profile.  Then,  a second  pass  is  started  which  uses  the  parameter  of  successfully  found  neighboring 
points  as  initial  values.  If  initial  conditions  are  found,  a Newton  iteration  is  performed  in  order  to  solve 
y = x(£,  T],  f) , which  uses  the  full  polynomial  approximation  of  x . 

3.2.3  Treatment  of  Wall  Points 


Points  at  walls  can  be  treated  differently  from  interior  ones,  since  a wall  point  in  the  CAA  grid  has  to  be  one  in  the 
RANS  grid  representation  of  the  wall  too.  Therefore,  the  search  for  initial  conditions  can  be  performed  only  over 
the  surfaces  of  the  RANS  blocks.  The  difficulty  which  arises  is  that  a wall  point  lies  at  a surface  of  a RANS  block 
and  the  Newton  iteration  for  the  determination  of  the  parameter  values  may  run  outside  the  block  which 
deterioriates  the  convergence  of  the  procedure.  Now,  however,  one  parameter  value  is  known  from  the  wall 
condition  and  consequently  the  search  can  be  performed  inside  the  face  of  the  block  (and  not  in  the  volume  of  the 
wall  cell).  In  this  case,  the  determining  vector  equation  y = x(£,T],£)  has  to  be  replaced  by  a set  of  two 
equations,  since  we  are  searching  two  parameter  values  only.  If,  for  example,  £ = const,  denotes  the  coordinate 
of  a wall,  ,rj  are  determined  such  that  the  difference  vector  y - x{f,7])  is  perpendicular  to  both  surface  tangent 

vectors  — and  — - — — . This  is  equivalent  to  the  condition  that  |y  - x(f,T])\  is  minimum.  This  approach 
dr] 

also  solves  the  problem  encountered  when  small  differences  between  the  wall  surfaces  of  both  grids  exist, 
originating  from  slightly  different  spline  approximations  in  both  cases. 
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3.2.4  Accuracy  of  the  Interpolation 

In  order  to  check  the  accuracy  of  the  interpolation,  the  test  function  f (x,  y)  = x2  + y 2 was  interpolated  from  the 

RANS  grid  onto  the  CAA  grid.  The  relative  error  (fCAA  - /*  )/ max(  fCAA  ) is  below  1 0-4  in  regions  where 
strong  mean  flow  gradients  are  expected,  cf.  Fig.4.  This  has  proved  to  be  sufficient  for  the  intended  purpose.  The 
error  displays  characteristic  extrema  where  large  jumps  in  gridspacing  on  the  RANS  grid  occur  or  the  grid  lines 
have  some  kinks  which  increase  the  error  in  the  spline  approximations.  It  should  however  be  noted  that  far  away 
from  the  body  the  RANS  solution  has  almost  no  gradients  anymore  and  in  the  vicinity  of  the  body  the  gridspacing 
is  very  small.  This  also  alleviates  accuracy  problems  of  the  interpolation. 


Fig.  4 The  relative  interpolation  error  on  the  CAA  grid 


4.  The  CAA  Calculations 

For  the  solution  of  the  linearized  Euler  equations  the  code  PIANO  was  used,  which  is  currently  under  develqpment 
at  DLR  [5,6].  The  LEE  solved  read  in  symbolical  flotation 


•f + v0-Vp  + vVp0+p0V-v  + pVv0=0 
at 

dv 

Po  ~^+  PoV0  •' Vv  + p0v  • Vv0  + pv0  • Vv0  = -Vp 
~ + v 0 • VP  + V •' Vp0  + TPo  v • V + jpV  ■ v0 = 0 


Here  p and  pare  the  density  and  pressure  disturbance  and  vis  the  vector  of  the  disturbance  velocity.  The 
quantities  with  index  0 are  the  mean  flow  values.  The  equations  were  brought  into  dimensionless  form  with  a 
reference  length  L , a reference  density  pw  and  a reference  velocity  which  was  choosen  to  be  the  velocity  of 
sound  far  away  from  the  body.  These  equations  are  solved  on  a structured  multi-block  grid,  like  the  one  used  in 
the  FLOWer  code.  The  spatial  derivatives  are  approximated  by  the  usual  4“’  order  DRP-stencils  [14]  and  time 
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integration  is  done  by  the  standard  4th  order  Runge-Kutta  procedure.  Inflow  and  outflow  boundary  conditions  were 
used  according  to  [14]. 

4.1  Filtering 

Since  the  utilized  DRP-Schemes  have  no  dissipation,  an  efficient  removal  of  spurious  waves,  i.e.  waves  which  can 
not  be  resolved  on  the  numerical  grid,  is  necessary.  The  usual  way  to  do  this  is  to  use  artificial  selective  damping 
[13].  This  approach  is  quite  expensive,  since  additional  terms  have  to  be  calculated  in  the  differential  equations.  A 
more  convenient  way  is  to  apply  periodically  during  the  calculation  suitable  digital  filters  on  the  fields  [12,15].  In 
PIANO,  symmetric  6Ul  and  8th  order  filters  are  implemented.  In  the  present  study,  the  filters  were  used  every 
hundred  time  steps.  The  6,h  ( N = 3)  and  8th  (N  = 4)  order  filter  arc  defined  as  [12] 

N 

Vi=a0 

j= i 

N = 3 : a0  = 0.6875  , o,  =0.46875  , a2  =-0.1875  , a,  =0.03125 
N = 4 : Oq  = 0.7265625  , a,  =0.4375  , a2  =-0.21875  , o3  =0.0625 
, a,  =-0.0078125 

where  \ffi  denotes  the  filtered  quantity  and  If/  the  unfiltered  one.  This  filter  is  applied  in  the  computational 
(£,  77,  C,)  -space  in  the  several  space  dimensions  subsequently.  No  filtering  was  used  in  the  first  N layers  normal 
to  a wall. 


4.2  Smoothing  Normal  to  Walls 

During  the  calculations  it  has  been  observed  that  very  strong  gradients  of  the  flow  variables  develop  in  the 
immediate  vicinity  of  walls,  especially  at  the  leading  edge  of  the  wing  section,  which  lead  to  numerical  instability. 
Examination  of  the  calculational  results  revealed  that  gridpoint  fluctuations  of  the  density  normal  to  the  wall  are  the 
driving  force  for  those  instabilities.  Therefore,  an  ad  hoc  procedure  was  applied  in  order  to  smooth  out  those 
gridpoint  fluctuations,  which  we  call  diagonal  smoothing'.  If  j denotes  the  direction  normal  to  the  wall  and 

j = 0 the  position  of  the  wall,  the  diagonal  smoothing  procedure  for  a flow  quantity  is  defined  as 

N 

¥'u=aoV'u+'La*<¥i-*j+*+V'M,i+t)  > J ~ 0,\,...,N  — ] 

k=\ 

The  coefficients  aj  are  the  same  as  for  the  filters  above.  The  value  of  N = 3 was  used  in  connection  with  6th 

order  filtering  and  N = 4 with  8th  order  filtering.  This  smoothing  procedure  was  not  applied  normal  to  the  trailing 
edge  face.  In  the  results  presented,  all  flow  quantities  were  smoothed,  although  in  some  calculations  it  would  have 
been  sufficient  to  smooth  the  density  alone. 

It  should  be  emphasized  however,  that  up  to  now  no  mathematical  justification  can  be  given  for  this  kind  of 
smoothing.  This  is  a problem  that  clearly  requests  further  investigations  in  the  future. 


43  Pressure  Pulse  Hitting  the  Leading  Edge 

A common  benchmark  problem  for  CAA  codes  is  the  scattering  of  a pressure  pulse  hitting  the  leading  edge  of  a 
wing  section.  Fig.5  shows  the  scattering  of  a pressure  pulse  located  at  x = -0.35,  y = 0 at  time  at  T = 0 in  front 
of  the  profile  for  the  RANS  mean  flow  and  8,h  order  filtering.  It  can  be  seen,  that  the  wave  front  travels  slightly 
faster  above  the  upper  side  of  the  profile  due  to  the  higher  velocities  of  the  mean  flow  there.  The  most  remarkable 
phenomenon,  which  was  not  observed  in  scattering  investigations  using  an  Euler  mean  flow,  is  however  the  strong 
generation  of  vorticity  in  the  boundary  layer  at  the  nose  of  the  profile.  This  vorticity  is  convected  by  the  mean  flow 
along  the  upper  surface  of  the  profile.  Fig.6  shows  the  contour  plot  of  the  v-velocity  component  and  the  local 
direction  of  the  velocity  field  by  means  of  momentary  streamlines  of  the  disturbance  velocity.  Fig.  7 shows  the 
corresponding  quantities  for  an  Euler  mean  flow  field,  where  almost  no  vorticity  is  generated  in  the  vicinity  of  the 
body. 
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c)  T=1.0  d)  T=1.5 

Fig.5  Scattering  of  a Pressure  Pulse  (RANS  mean  flow) 
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Fig.6  T=1.5,  RANS  mean  flow:  v-velocity  field  with  streamlines  of  the  disturbance  velocity 
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Fig.  7 T=1.5,  Euler  mean  flow:  v-velocity  field  with  streamlines  of  the  disturbance  velocity 

A physical  explanation  for  this  strong  generation  of  disturbance  vorticity  in  regions  where  the  mean  flow  field  has 
vorticity  maxima'  too  can  be  given  if  one  considers  the  compressible  vorticity  transport  equation: 

Dm  „ . V/txVo 

= -o)Vv  + (w- V)vh — ---  - 

Dt  p- 

After  linearization  the  first  term  on  the  right  hand  side  gives  a term  consisting  of  the  product  of  the  vorticity  of  the 
mean  flow  and  the  divergence  of  the  disturbance  velocity.  Since  the  mean  flow  vorticity  takes  huge  values  in  the 
boundary  layer  this  results  in  a strong  source  term  for  the  disturbance  vorticity  if  the  pressure  wave  arrives.  Similar 
vortex  structures  have  been  found  by  Manoha  et  al.  [10]  performing  a large  eddy  simulation  of  the  flowfield 
around  a NACA  0012  airfoil. 


Fig.  8 T=1.5,  Pressure  on  a circle  around  the  profile  for  different  mean  flow  fields  and  filters 
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In  order  to  compare  the  influence  of  the  different  mean  flow  fields  and  two  different  filters,  for  T=1.5  the  pressure 
on  a circle  (origin  (0.5, 0.0),  radius  0.8)  around  the  profile  was  plotted,  cf.  Fig.8.  The  angle  0°  corresponds  to  the 
downstream  direction. 

It  can  be  seen  that  the  difference  between  the  6th  and  the  8th  order  filters  is  much  larger  than  the  difference  between 
RANS-  and  Euler  mean  flow.  As  one  would  expect,  one  recognizes  that  the  pressure  profile  for  the  8th  order  filter 
has  sharper  gradients  than  that  of  the  6th  order  filter.  The  largest  difference  between  the  two  filters  occurs  near  the 
trailing  edge  of  the  profile  (domain  near  angle  0°),  which  becomes  clear  if  on  considers  that  there  the  sharpest 
gradients  in  the  disturbance  flow  field  develop.  The  very  small  difference  between  the  pressures  of  the  RANS  and 
Euler  mean  flow  fields  indicates  that  almost  no  sound  is  generated  by  the  vortex  structures  mentioned  above  and 
that  the  influence  of  mean  flow  field  wall  boundaty  layers  can  be  neclected  in  the  scattering  problem  considered. 


5.  Conclusions 

For  the  solution  of  the  linearized  Euler  equations  with  a RANS  mean  flow  field  the  RANS  grid  can  normally  not 
be  used.  Therefore  an  interpolation  procedure  for  the  mean  flow  field  between  curvilinear  structured  multi-block 
RANS-  and  CAA-grids  was  suggested,  which  is  based  on  a higher  order  polynomial  approximation  of  the  grid 
variables  on  each  cell.  The  coefficients  of  the  polynomial  are  determined  from  the  function  values  and  its  2nd 
derivatives  at  the  comers  of  the  cell.  The  2nd  derivatives  are  determined  by  one  dimensional  cubic  spline 
interpolation  along  the  gridlines.  Of  crucial  importance  is  the  provision  of  suitable  initial  conditions  for  the  iterative 
determination  of  the  the  parameter  values  of  the  grid  functions. 

The  interpolation  procedure  was  applied  on  the  scattering  problem  of  a pressure  pulse  hitting  the  leading  edge  of  a 
wing  section  at  an  angle  of  attack  of  CC  — 1°.  Remarkable  is  the  strong  production  of  vorticity  in  the  nose  region  of 
the  profile.  This  vorticity  is  transported  by  the  mean  flow  along  the  upper  surface  of  the  profile,  but  obviously  does 
not  produce  sound  waves,  as  shows  the  comparison  with  calculations  using  an  Euler  mean  flow.  In  the  scattering 
problem  considered,  the  differences  in  the  pressure  fields  of  the  RANS  and  Euler  mean  flow  cases  are  very  small. 
The  difference  between  6th  and  8th  order  filtering  are  most  pronounced  in  the  wake  of  the  profile,  where  sharp 
gradients  evolve  which  were  more  intensely  smoothed  by  the  6lh  order  filter. 

An  open  problem  that  remains  is  the  strong  mathematical  justification  of  the  diagonal  smoothing  procedure,  which 
had  to  be  applied  in  immediate  vicinity  of  solid  walls,  in  order  to  stabilize  the  calculations. 
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Question: 

The  RANS  solution  is  converged  in  the  sense  that  mass  and  momentum  balances  on  the  RAND 
cells.  Do  you  check  that  after  interpolation  the  interpolated  mean  flow  satisfies  mass  and  momentum 
balance  on  the  CAA  grid?  If  not,  the  linearized  Euler  equations  you  solve  in  a discretized  sense  may 
contain  artificial  diffusive  fluxes  other  than  zero  on  their  right  hand  sides. 

Answer: 

The  balance  of  the  interpolated  fluxes  is  not  checked.  Since  the  linearized  Euler  equations  are 
solved  in  nonconservative  form,  artificial  fluxes  are  always  present  on  the  right  hand  side:  but  are 
assumed  to  be  small.  We  assume  that  the  modeling  error  in  the  RANS  calculations  is  at  least  one  order  of 
magnitude  larger  than  possible  errors  introduced  by  the  interpolation. 


